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ABSTRACT 

We assess the relative merits of weak lensing surveys, using overlapping imaging data from the 
ground-based Subaru telescope and the Hubble Space Telescope (HST). Our tests complement similar 
studies undertaken with simulated data. From observations of 230,000 matched objects in the 2 square 
degree COSMOS field, we identify the limit at which faint galaxy shapes can be reliably measured 
from the ground. Our ground-based shear catalog achieves sub-percent calibration bias compared to 
high resolution space-based data, for galaxies brighter than i 1 ~24.5 and with half-light radii larger 
than 1.8". This selection corresponds to a surface density of 15 galaxies arcmin -2 compared to 
~ 71 arcmin -2 from space. On the other hand the survey speed of current ground-based facilities is 
much faster than that of HST, although this gain is mitigated by the increased depth of space-based 
imaging desirable for tomographic (3D) analyses. As an independent experiment, we also reconstruct 
the projected mass distribution in the COSMOS field using both data sets, and compare the derived 
cluster catalogs with those from X-ray observations. The ground-based catalog achieves a reasonable 
degree of completeness, with minimal contamination and no detected bias, for massive clusters at 
rcdshifts 0.2 < z < 0.5. The space-based data provide improved precision and a greater sensitivity to 
clusters of lower mass or at higher redshift. 

Subject headings: cosmology: observations - gravitational lensing - instrumentation 



1. INTRODUCTION 

Dark matter dominates the gravitational component of 
the cosmic energy density and thus provides the frame- 
work for structure formation in the Universe. How- 
ever, by its nature, the distribution and cosmic growth 
are challenging to observe. The most promising probe 
is weak gravitational lensing: analysis of the distorted 
shapes of ordinary galaxies behind foreground mass con- 
centrations. Several numerical techniques are now avail- 
able to recover the projected mass distribution from 
these distortions, and tests on sim ulated datasets are un- 
derway to verify th eir precision (jHevmans et al. 2 005a: 
iMassev et al. 2007d ). There is great optimism in the 
weak lensing community that such methods will enable 
both the tomographic mapping of dark matter structures 
in time and space. This will also provide a robust statis- 
tical mea sure of the nature of dark ene rgy over redshifts 
0< z <1 (jMellier 1999t iRefregier 20031 ) 

Observational progress has been particularly 
dramatic. The first detections of statistical 

"cosmic shear" were only pu blished in 2000 
dBacon. Refreeier fc Ellis 200a _ iKaiser et al. 20001 : 
IWittman et al. 20001 : Ivan Waerbeke et al. 200(1 . 

In the subsequent seven years, weak lensing sur- 
veys ha ve measured the dark matter power spec- 
trum dBrown et al 20031: iHevmans et al. 2005bt 
iHoekstra et al. 20061 Scmbol oni et al. 2006f). traced 
the evolution o f structure dBacon et al. 20051 : 
iKitching et al. 20071 : IMassev et al. 2007al ). enabled 
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the co nstruction of lensin g-sel ected cluster cata- 
logs dMivazaki et al. 2002al: IWittman et al. 20061 : 
ISchirmer et al. 2007t iMivazaki et al. 20071 ). and 
non-parametrically reconstructed the total mass 
distribution bot h in clusters (|Kneib et al. 20031 
IClowe et al. 20061: LTee et al 2007al ) and on larger 
scales (jMassev et al. 200 7b). As a result, weak lens- 
ing has been identified as the most promising route 
to understanding the nature of dark energy by the 
ESA-ESO Working Group on Fundamental Cosmology 5 , 
joint NSF-NASA-DOE Astronomy and Astrophysics 
Advisory Committee 6 , and NSF-DOE High Energy 
Physics Advisory Panel Dark Energy Task Force 7 . 

The primary signal of any weak lensing analysis is the 
statistically coherent distortion of background galaxies 
along adjacent lines of sight. The main sources of statis- 
tical noise are the finite density of galaxies that can be 
sufficiently well-detected and resolved for accurate shape 
measurement, plus their intrinsic morphologies. The 
density of resolved galaxies also governs the angular res- 
olution and fidelity of a reconstructed mass map which, 
in turn, determines the limiting halo mass that can be 
detected. On the other hand, statistical analyses of the 
dark matter power spectrum are less concerned with in- 
dividual halos but require panoramic fields to counter 
the effects of cosmic (sample) variance. Minimizing sta- 
tistical errors in such an analysis, within a finite survey 
lifetime, requires an optimal balance between area and 
depth. 

A key debate in the development of future weak lensing 
experiments concerns the relative merits of ground- ver- 
sus space-based platforms. Ambitious surveys now being 

5 http: //www. stecf . org/coordination/esa_eso/cosmology .php 

6 http : //www . nsf . gov/mps/ast/aaac . j sp 

7 http : //www . nsf . gov/mps/ast/detf . j sp 
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planned with dedicated, ground-based facilities (eg VST- 
KIDS, DES, Pan-STARRS, LSST). These are driven by 
technological progress including panoramic cameras with 
small optical distortions, highly sensitive imaging detec- 
tors, and (in the case of Pan-STARRS) on-chip active 
correction to reduce the width of the point spread func- 
tion (PSF). Future surveys spanning significant fractions 
of the celestial sphere are envisaged, promising tight con- 
straints on the cosmological parameters. 

However, measurements with current ground-based fa- 
cilities are limited by the size and temporal variations of 
the PSF. There is concern in many quarters that wide- 
field facilities operating in space (e.g. DUNE, SNAP, 
JDEM) will ultimately be required to achieve the pre- 
cision required (particularly) to distinguish between var- 
ious models of dark energy. Space-based facilities will 
be more costly but will likely offer increased depth, bet- 
ter photometric performance and a stable PSF. The key 
issue in gauging their merits is not statistical error, but 
the extent to which potential biases in ground-based data 
may act as a "systematic floor" to prevent complete ex- 
ploitation. 

Some valuable answers can be obtained by com- 
paring simulat e d ground and spa ce-based images, 
i)Wittman 2005t lLampton et al. 20061) and the Shear 
TEsting Programm e (STEP: iHevmans et al. 2005al : 
iMassey et al. 2007cl ). However, the input parameters 
used to generate the simulated data may not be real- 
istic or address all the instrumental idiosyncrasies. Of 
particular concern are the stability and vagaries of the 
PSF. No simulations have yet adequately addressed this 
point - which may, ultimately, be the limiting problem 
for ground-based data. It is often argued that future fa- 
cilities will be carefully designed to mitigate any limita- 
tions realized with current observational facilities. While 
progress can no doubt be expected, both on the ground 
and in space, we believe many lessons can be learned 
from extant data and hardware with proven engineering 
pedigree. 

In this paper, we present the first direct comparison 
of weak lensing analysis for the same sky field using 
ground and space-based data. Deep, panoramic imag- 
ing has been obtaine d for the 1.64 deg 2 COSMOS field 
(jScoville et al. 200~7at ) by both the Advanced Camera for 
Surveys (ACS) on board th e Hubble Space Telescope 
(HST) (|Scoville et al. 2007H) and the Suprime-Cam im- 
ager at the prime focu s of the Subaru 8.2m telescope 
(jTaniguchi et al. "20 07a). In both cases, the entire field 
was covered by mosaicing many independent exposures. 
The SuPrimeCam instrument was constructed with weak 
lensing analysis particularly in mind, and currently pro- 
vides the best image performance available from any 
ground-based telescope, in terms of optical distortions 
over a large field. A comparison of these datasets should 
therefore provide a realistic and valuable assessment of 
the relative performance of state-of-the-art imagers on 
the ground and in space. 

The paper is organized as follows. In Sj2j we briefly re- 
view the relevant theory. In we describe the two data 
sets, data reduction pipelines and weak lensing analyses. 
We then present the results. In 21 we compare shear 
measures on a galaxy-by-galaxy basis to determine the 
optimum depth at which the ground-based data matches 
the performance of the (deeper) space based data. This 



permits us to determine the relative survey speeds of 
Subaru and HST for high precision cosmic shear exper- 
iments. In fJSl we construct maps of the mass distribu- 
tion, treating the Subaru and HST maps as independent 
probes of the same field, and contrast these against X-ray 
data. This permits us to evaluate the completeness and 
reliability of a lensing-selected halo catalog, and evalu- 
ate the precision of their inferred masses as a function 
of redshift. In we summarize our results and discuss 
their wider implications for future missions. 

2. REVIEW OF WEAK LENSING THEORY 

Gravitational lensing by foreground mass structures 
distorts an image plane of distant galaxies /(x) via a 
coordinate transformation 
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where 5xi(x) is the deflection angle of the light rays. The 

convergence 



k(x) 



2c 2 



g(z)p(>i 1 z)d.z , 



(2) 



describes overall dilations and contractions. It is propor- 
tional to the total mass density p projected along a line 
of sight, where the lensing sensitivity function 
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reflects the efficiency of foreground gravitational lenses 
at different redshifts - containing a ratio of the angular 
diameter distance to a lens, the background source, and 
between the two. This can be more simply written as 
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in terms of a 2D, projected version \P(x) of the Newto- 
nian gravitational potential. Two components of shear 
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describe stretches and compressions along (at 45° from) 
the x-axis. 

The observed shapes of background galaxies can be 
described by combination of their Gaussian-weighted 
quadrupole moments 
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Although k is generally the desired quantity, and could 
be obtained in principle from measurements of galaxy 
sizes © or fluxes, this has proved difficult in practice, 
because expectations for these quantities prior to lensing 
are unknown. On the other hand, while galaxies have 
a natural dispersion of intrinsic ellipticitics |(7J) , they are 
(almost) uncorrelated with each other in the absence of 
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lensing, i.e. (e^) = 0. Any correlation between the ob- 
served ellipticities of galaxies seen along adjacent lines of 
sight arises because their light has traversed similar in- 
tervening large-scale structure p(x, z). In practice, cor- 
rections to measured ellipticities also need to be made 
for the smearing of galaxies by the PSF, and for the dif- 
fering susceptibilities of some galaxy morphologies to an 
input shear. For more details of thi s procedure, see e.g. 
(|Kaiser. Squires fc Broadhurst 19951 ). 

The observed shear can finally be transformed into con- 
vergence through their close relation in Fourier space 

. (if - t&n + (2Wj| m 
WTW) (9) 

([Kaiser fc Squires"T 993). This is typically some amount, 
to reduce noise. Furthermore, like any scalar quantity ex- 
tracted from a vector field, a convergence signal can also 
be split into two indepen dent components, n — k e + i k b 
(jKing fc Schneider 200 J ). The grad-like "E-mode" is the 
signal produced by weak lensing. The curl-like "B-mode" 
is not produced by ph ysical processes (except at very low 
levels, as described bv lSchneider et al. 2002| ). and there- 
fore ought to be consistent with zero in the absence of 
systematics. Usefully, it contains the same noise proper- 
ties as the -E-mode signal - so it acts as an independent 
realization of noise in the field, and any significant devi- 
ations from zero alert to the presence of residual system- 
atics (such as imperfect correction for the PSF). 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. The COSMOS Data Sets 

Our data all cover the COSMOS survey field, a 
1.64 deg 2 contiguo us square, centered at 10:00:28.6, 
+02:12:21.0 (J2000) (|Scoville et al. 2007al ). The ground- 
based imaging was obtained in eleven mosaiced pointings 
of the Suprime-Cam camera at the prime focus of th e 
Subaru telescope on Mauna Kea ( Mivaza ki et al. 20 02b). 
These were taken on the 18th and 21st of February 2004, 
nights selected for their excellent observing conditions: 
the mean seeing was 0.54" ± 0.03". The field consti- 
tutes part of a larger weak lensing survey discussed, along 
wit h full details of the p rimary data reduction pipeline, 
in (|Mivazaki et al. 20071 Green et al. in prep.). In fact, 
the relevant field in that survey covered a slightly larger 
area than the COSMOS field. The Subaru imaging was 
truncated when matching galaxy catalogs and was trun- 
cated after making convergence maps, to avoid edge ef- 
fects associated with the Fourier transform operations in 
equation ©. 

Our comparison is made possible by the unique avail- 
ability of deep, panoramic space-ba sed imaging of the 
COSMOS field (|Scoville et al. 2007H) . During HST cy- 
cles 13 and 14, 577 slightly overlapping pointings were 
obtained from the Advanced Camera for Surveys (ACS) 
on board the Hubble Space Telescope. Four dithered ex- 
posures at each pointing were stac ked using the drizzle 
algorithm (jFruchter fc Hook 20021) to improve the native 
pixel scale of 0.05" and recover a final pixel scale of 0.03". 
Full details of the primary data reduction pipeline for the 
HST images are given in (|Ko ekemoe r et al. 20071 ) . 

It is important to emphasize that both the ACS and 
Suprime-Cam data exhibit idiosyncrasies that present 
significant challenges for weak lensing analysis. For ex- 
ample, the atmospheric seeing varied during the two 



nights over which the Subaru data were obtained; and 
the distortions of the telescope's primary mirror under 
a gravity load were only passively corrected via a look- 
up table as it followed the field. In a future ground- 
based experiment, such as LSST or Pan-STARRS, see- 
ing variations could be normalized over a survey by 
stacking a very large number of short, independent ex- 
posures taken over a long time period. Dome seeing 
could likewise be improved with future technologies. 
And while the telescope superstructure is particularly 
rigid at Subaru, active correction of the mirror sup- 
port could undoubtedly improve future designs. Equiv- 
alently, the sky background seen from HST is affected 
by Earthshine that de pends on the telescope pointing 
(jLeauthaud et al. 2007fh The HST PSF also varies over 
time due to thermal fluctuations during ea ch low-Earth 
orbit pee et al. 2007at iRhodes et al. 20071 ). Finally, the 
charge transfer efficiency of the ACS CCD detectors had 
been significantly degraded by high energy particles by 
the time the COSMOS data w ere obtained, and w orsened 
during the observing window (|Rhodes et al. 20071 ). None 
of these problems are inherent to all space-based obser- 
vations: future missions might minimize or eliminate all 
three effects by adopting a regular observing pattern, 
orbiting the Lagrange point L2, and using radiation- 
hardened CCDs. However, the weak lensing analysis of 
existing space-based data is indeed compromised by the 
extent to which such hardware variations can be mod- 
eled. In this sense, our comparison is actually more infor- 
mative than one based on simulated data that reproduces 
only idealized and mean instrumental characteristics. 

The relevant characteristics of the two data sets are 
summarized in Table [IJ including limiting depths for a 
point source at 5cr, in a 3" ape rture from the gro und 
and a 0.15" aperture from space (|Capak et al. 2007[) . In 
addition to these images, the COSMOS field has been 
observed across all wavelengths from radio to X-rays. 
Of particular relev ance here are (i) de ep X-ray obser- 
vations by XMM (Hasingc r et al. 2007t ). which can be 
used to locate massive structures via thermal emission 
from hot gas; and (ii) multicolor optical and near-IR 
imaging campaigns from the Subaru, Canada France 
Hawaii, Cerro Tololo and Kitt Peak telescopes, which 
provide 15 additio n al bands and photom etric redshifts 
ijCapak et al. 20071 : iMobasher et al. 20071) . The pho- 
tometric redshift estimation code uses Bayesian priors 
based on an adopted luminosity function, and includes 
reddening based on both Galactic and Calzetti extinc- 
tion laws. The results were calibrated using 868 galaxies 
in the field brighter than i' — 24 and with spectroscopic 
redshifts. For galaxies closer than z = 1.2, the rms scat- 
ter in (z phot - z spec ) I (1 + z spec ) is 0.031. 

3.2. Object Detection 

Objects were detected in the Subaru images using 
hfindpeaks from the IMC AT package 8 . This finds the 
centroid and scale size r g that maximizes the peak S/N 
of the image after smoothing with a Gaussian. The code 
also returns the half-light radius, r^, of each galaxy. 
Galaxies were initially detected to magnitudes fainter 
than those for which it is possible to accurately measure 

8 Nick Kasiser's imcat software package is available from 
http : //www. if a.hawaii . edu/~kaiser/imcat/ 
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TABLE l 

Survey Characteristics 





\ -1 v*i^i i vi r\ 

Vjl U L111L1 


Spucc 


Instrument 


Subaru / Suprime- Cam 


HST/ ACS 


Primary aperture 


8.2m 


2.4m 


Exposure time 


5 x 360s 


4 x 507s 


Total survey time 


5 hours 


325 hours 


Filter 


i 1 


F814W 


Limiting AB magnitude 


26.2 


26.6 


Field of View 


2.14 deg 2 


1.67 deg 2 


Pixel Scale 


0.202" 


0.03" 


Point Spread Function 


0.68" 


0.12" 



shapes. To reduce noise in the final analysis, weights 
were given to each galaxy, and galaxies with a detection 
S/N< 14 were removed from the catalog altogether. The 
resulting surface density is n ga i = 42 galaxies arcmin~ 2 . 

Objec ts were detected in the ACS images using SEx- 
tractor (iBertin fc Arnouts 19961) in a dual Hot/Cold con- 
figuration (jLeauthaud et al. 20071 ) . designed to identify 
both large and small objects while avoiding fragmenta- 
tion of the former, or merging of the latter. The SExtrac- 
tor centroids were then improved, and the best-fitting 
scale size was selected, via an iterative process during 
shape measurement. Galaxies smaller than d = 0.11" or 
fainter than S/N « 20 were removed from the catalog, 
and a weighting scheme was app lied to faint galaxies as a 
function of their detection S/N (jLeauthaud et al. 2 007) . 
Note that an absolute calibration of the S/N was diffi- 
cult to determine in practice, because flux in adjacent 
pixels becomes correlated during DRIZZLE. The S/N 
cut corresponds approximately to a limiting magnitude 
F814IF(AB) < 26.5 for a point source. The resulting 
surface density is n ga i — 71 galaxies arcmin , with a 
median redshift z mc( j = 1.2. 

3.3. Shear Measurement 

Because the image characteristics of the two data 
sets are quite different, we adopted separate methods 
to measure galaxy shapes, remove PSF effects, and 
ultimately obtain the weak lensing shear signal. Each 
of these methods has been optimized for the respective 
data sets, so our comparison will necessarily incorporate 
the limitations of each pipeline. We believe this is in 
the spirit of a fair comparison of ground versus space. 
To minimize any differences arising entirely from the 
algorithms themselves however, we have intentionally 
adopted related methods from the same generation 
of software development and codes that have been 
well- tested. Alth o ugh newer shear measurement meth- 
ods dKaiser 20001: iBridle et al. 20011 iDahle et al. 2002 : 



IBernstein fc Jarvis 2002 
" TKuiiken 2006 



Refregier & Bacon 200 
Hira ta fc Seliak 200, 

Massev et al. 2007dt iNakajima k, Bernstein 2007 ) 

may offer improved performance, none has yet been 
sufficiently tested across both observing regimes. 

The Subaru im ages were analyzed with the 
(lHamana et al. 2003T) implementation of the widely 
used (|Kaiser. Squires k, Broadhurst 19951 hereafter 
KSB) shear measurement method. This partic- 
ular implementation is a derivative of the "LV" 
pipeline tested in the Shear TEstin g Program 
(jHevmans et al. 2005at IMassev et al. 2007cD . 



The HST ima ges were analyzed with the 
([Rhodes et al. 200"lL hereafter RRG) shear mea- 
surement method. This is a perturbation of the KSB 
method for space-based data. It calculates the same 
quadrupole moments, but corrects them individually 
for the effects of convolution with the PSF, and only 
in the final stage takes the ratio ([7]). This is neces- 
sary because the small and cuspy diffraction-limited 
PSFs otherwise introduce divisions by very small (and 
noisy) numbers. R RG been applied to HST W FPC2 
dRhodes et al. 200 ID. STIS dRhodes et al. 20041) . and 
A CS data (jMassev et al. 2007af l. The ACS pipeline was 



thoroughly tested on simulated images during the cre - 
ation of the COSMOS catalog (jLeauthaud et al. 20071) . 
and also for a continuation of STEP using simulated 
space-based images. 

4. STATISTICAL APPLICATIONS 
4.1. Shear-shear comparisons 

We shall now compare the global properties of our 
ground- and space-based shear catalogs, to determine the 
depth (and galaxy surface density) at which reliable shear 
measurement is possible from the ground. This will be 
relevant for many statistical applications, including mea- 
surements of the angular shear-shear correlation function 
that are typically used to constrain cosmological param- 
eters. In such analyses, where statistical noise is reduced 
by averaging over many lines of sight, the key issue is the 
reliability and level of residual systematics in the shear 
measurement. 

We asses the performance of the ground-based shear 
measurements against those from the same galaxies in 
space-based data, making the necessary but reasonable 
assumption that the shapes are much more reliable when 
measured from the much higher resolution images with 
a smaller PSF. Such a comparison is clearly only pos- 
sible for the subset of galaxies contained in both cata- 
logs. The two quantities of interest will be linearity in 
the comparison (the slope of the shear-shear comparison 
is equivalent to the STEP "calibration bias" parameter 
m) and the scatter (which represents the combined shear 
measurement noise from both HST and Subaru, plus any 
systematic effects). 

We match galaxies whose positions agree to within 1", 
and produce a common catalog containing n ga i = 32 
galaxies arcmin -2 . Many objects in the Subaru galaxy 
catalog without matched counterparts in HST galaxy 
catalog have half-light radii on the limits of seeing and are 
likely to have been revealed as stars by the higher reso- 
lution data; in any case, the omitted galaxies had below- 
average weights in the Subaru catalog. The remaining 
unmatched objects are a combination of noisy/skewed 
galaxies with offset centroids, or galaxies that lie in re- 
gions of the HST images masked because of scattered 
light from nearby bright stars. For the following tests 
tests, we shall ignore the weights on remaining galaxies, 
and treat all objects equally. 

Figure [1] shows a comparison of the shear signal for 
the matched galaxies. Since errors are present in both 
axes, to calculate the best-fit linear relationship we adopt 
a least squares method that minimizes the perpendicu- 
lar distance to the best fit line (instead of one that as- 
sumes one variable is error-free). Ideal shear measure- 
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Fig. 1. — Comparison of shear measured from galaxies seen in 
both Subaru and Hubble Space Telescope images. The greyscale 
shows the number of galaxies with different shear measurements. 
The outer contour includes 90% of the galaxies, and successive in- 
ner ones include 10% fewer. The solid line is the least squares linear 
relation. Its slope of 0.87 indicates that shears have been under- 
estimated in the ground-based analysis, or that the catalog is still 
partially contaminated by stellar sources. This value is insensitive 
within 0.01 to the reintroduction of galaxy weights. Furthermore, 
the non-Gaussian wings of the scatter extend well beyond the rms 
error of 0.16, shown as dashed lines. 

ment from both instruments would yield a best-fit slope 
of unity. There will inevitably be a small amount of scat- 
ter, because the weight function ([8]) is not necessarily 
the same size r g in the ground- and space-based analy- 
ses. In practice, we find a best-fit slope of 0.87, indicat- 
ing that the shears have been underestimated from the 
ground. The measurement noise is also problematic, with 
er 7 = 0.16 per component (perpendicular to the best-fit 
line; note that this does not include intrinsic source ellip- 
ticity variance because the galaxies are matched) and a 
skewed, non-Gaussian distribution of outlying shear es- 
timates that would render a cosmic shear analysis less 
stable. 

The overall performance in figure[T]is a superposition of 
good shears from bright and (in particular) large galax- 
ies, plus smaller objects that cause most of the bias and 
scatter. Indeed, systematic errors could be completely 
eliminated by using only the very largest galaxies. How- 
ever, the statistical noise in a cosmic shear analysis of 
shear-shear correlation functions scales as cr 7 /^/7i ga i. An 
optimal strategy for any particular ground-based survey 
will involve catalog cuts requiring a trade-off between 
systematic and statistical errors. However, the optimal 
cuts will vary as a function of survey area and depth. To 
produce a result of general interest, we therefore show 
in Figures [2] and [3l the resulting calibration bias, scatter 
and galaxy density for a range of possible cuts in galaxy 
size and magnitude. 

A simple result emerges from Figure O It is notice- 
able from the horizontal and vertical contours that size 
and magnitude cuts seem to neatly parametrize inde- 
pendent sources of error. Using existing shape measure- 
ment methodology, shear can be measured from galaxies 
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Fig. 2. — Relative calibration between shear measurements from 
galaxies in Subaru and Hubble Space Telescope data, for galaxies of 
different sizes and magnitudes. The contours show deviations from 
a slope of unity in figures [T] and [4] which would have been ideal. For 
faint galaxies, these deviations tended to be an underestimation in 
the Subaru pipeline relative to HST. There is some evidence that 
shears are overestimated in large, bright galaxies, although the 
small number of these objects means that the extrapolation is less 
certain. The calibration biases are calculated locally, for galaxies 
only in a given cell of {size, magnitude} space. On the other hand, 
the grey-scale shows the cumulative number density of galaxies 
n ga [ that would remain in a ground-based catalog, were cuts to 
be applied at the local values (i.e. including all larger and brighter 
galaxies). 
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Fig. 3. — Combined noise from shear measurements of galax- 
ies matched in catalogs from Subaru and Hubble Space Telescope 
data. The contours show cr 7 as a function of galaxy size and mag- 
nitude. As in figure [2] these are calculated only for galaxies with 
that particular size and magnitude. The contours close at the top 
merely because there are very few large, faint galaxies, so the rms 
scatter increases. The grey scale again shows the total number 
density of available galaxies. 
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Fig. 4. — As for figure [T] but for the subset of galaxies brighter 
than i' = 24.5 and larger than = 1.8". The total least squares 
slope is 0.97, implying an almost unbiased recovery of the shear 
signal from Subaru, and the data is better and more symmetrically 
enclosed within the rms scatter of 0.11. 

brighter than i' = 24.5 and larger than — 1.8, with 
measurement noise cr 7 ~ 0.03 and a calibration bias less 
than 3% (and only 1% with galaxy weighting), which is 
accep table for competitive constraints from future sur- 
veys (Rcfrcgic r et al. 20041 ) . This leaves a surface den- 
sity of n ga i = 15 galaxies arcmin 2 from the ground, with 
a median redshift of z mc d =0.8. The comparison with 
space-based data for these cuts is shown in figured] 

Note that we have not been able to test the reliability 
of space-based shear measurements using this method, 
nor even considered the population of small galaxies re- 
solved only from space. Without data even better than 
ACS imaging to compare to, we resort to simulations. 
The full R RG pipeline was calibra ted against simulated 
images by (|Leauthaud et al. 20071 ). However, of the 71 
galaxies arcmin -2 successfully in r eal A CS images, only 
the brightest 40 could be used by (jMassev et al. 2007af ) 
to minimize the B-mode signal and overcome problems 
of CCD charge transfer inefficiency (CTI). This limita- 
tion clearly needs to be overcome: perhaps via a CTI 
correction algorithm like that developed for STIS by 
(|Bristow et al. 20041 ). and radiation- hardened detectors 
in future telescopes. 

4.2. Survey speed 

Our ground versus space shear comparisons have im- 
portant implications when it comes to considering the 
optimal approach for measuring shear for cosmological 
applications. Although the etendue of instruments can 
be expected to increase both on the ground and in space, 
we will base our discussion on the imaging depth and 
fields of view of our HST and Subaru surveys. 

An important criterion is what can be accomplished 
in a given amount of observing time; the HST and Sub- 
aru requirements for our comparison are summarized in 
Table [lj . HST overheads approximately tripled the on- 
source exposure time, and, at Subaru, high quality imag- 
ing was secured during what might be considered a fortu- 



itous observing window. Noting that (|Bacon et al. 2001h 
found images with seeing worse than ~ 0.8" of little use 
for weak lensing analysis, coupled with observational vis- 
ibility, it seems reasonable to incorporate a factor of at 
least four inefficiency for a generic survey: even for a 
superb ground-based facility such as Subaru, on an ex- 
cellent site such as Mauna Kea 9 . Based on its superior 
field of view, Subaru is then ~ 24 times faster than HST 
in useful mapping speed. As, to first order, the signal 
to noise in statistical analyses increases as ^Jn ga \, for a 
fixed survey lifetime, this corresponds to a ~ 5-fold im- 
provement in signal to noise. 

This simplistic analysis is of course mitigated by the 
higher resolution available from space. We next insert 
the gain in surface density, viz 71 galaxies arcmin -2 
resolved in our space-based imaging c.f. 15 arcmin -2 
from the ground. We will assume that the additional, 
small galaxies have a similar distr ibution of intrinsic 
ellipticities as the lar ger ones (c.f. iMassev et al. 2004 
iLeauthaud et al. 2007h and that the measurement noise 
on an average survey galaxy is constant (since the size 
distribution of resolved galaxies compared to the PSF 
size is roughly independent of the PSF size). Incor- 
porating this increased background surface density, the 
ground-based gain per unit time drops to only a factor 
of 2.3. 

Equally important to the increased density of galaxies 
are their higher rcdshifts. Distant galaxies are more sen- 
sitive to low-redshift lenses, and sensitive to more total 
lenses. The shear signal grow s proportionally to the me- 
dian source redshift as (I.Iain fc Seliak 1997h . With 
the redshift distributions for galaxies shown in figure [5] 
the total gain in signal to noise for a 2D weak lensing 
survey conducted from Subaru over one conducted with 
HST is only about 1.7. 

Perhaps the most important advantage of space is the 
increased redshift range of resolved galaxies. This better 
enables their stratification into redshift bins for tomo- 
graphic (3D) analyses. Deep infrared imaging, needed 
for accurate photometric redshifts, is also likely to re- 
main the province of space-based observatories. To- 
mographic techniques can tighten the constraints on 
cosmological p arameters Am and <jg by a factor of 
at least three (IMassev ct al . 2007al) and potent ially as 
much as five ([Heavens. Kitching &: Tavlor 20061 ). Fur- 
ther advantages of these techniques includes the elim- 
ination of unwanted signal f rom adjacent galaxies ' 
intrinsically-correlated sha pes ([King fc Schneider 20031 : 
iHevmans fc Heavens 20031 ) . While wide-field ground- 
based instruments may therefore yield significant im- 
provements for Dark Energy Task Force "Stage 3" sur- 
veys, advanced analysis techniques for "Stage 4" surveys 
will realistically be possible only with space-based facili- 
ties. These will bring new scientific opportunities, cross- 
checks for systematics, and greater efficiency. 

5. MASS MAPS AND HALO DETECTION 



9 http: //www. cf ht . hawaii . edu/Instruments/Imaging/Megacam/ 
observingstats.html. future surveys such as Pan-STAKKS and 
LSST, which plan to co-add many short exposures with indepen- 
dent PSFs, may achieve near-uniform image quality by rejectng a 
certain fraction of exposures. But the relevant figure of merit is 
still the fraction of time spent with seeing better than 0.8". 
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Fig. 5. — The dotted lines show the redshift distribution of 
source galaxies from Hubble Space Telescope (black) and Subaru 
(red) imaging surveys, for the catalog cuts used in figure [4] The 
jagged lines show the measured photo-zs, and the smooth curves 
assume a simple para metric form for the backgroun d galaxy red- 
shift distribution from (Smail. Ellis & Fitchett 1994), with a = 2, 
j3 = 1.5, 2 me d =0.8 or 1.2, and an overall normalization to repro- 
duce the observed number density of galaxies. The solid lines show 
the corresponding lensing sensitivity functions calculated from the 
analytic curves. Thees lie always in front of the source galaxies 
but are notably higher for a space-based survey, particularly at 
redshifts greater than 0.5. Extending this redshift coverage is cru- 
cial for 3D tomographic analysis techniques to measure the growth 
of structure. 

We now investigate the reconstruction of maps of the 
mass distribution (figures [6] and [7]), and the detection 
of individual mass peaks. The mass and redshift dis- 
tribution N(M, z) of several thousand lensing-selected 
clusters could be used to constrain cos mological models 
(|Hamana et al. 2003t IWang et al 2004 Green et al. in 
prep.). Additionally, the physical properties of the dark 
matter particles can be investigated by comparing the 
detailed distribution of dark matter with that of baryons 
(|Clowe et al. 20061: 1 Jee et al. 2007bD . The key issues will 
be the angular resolution of reconstructed mass maps, 
as well as the mass and redshift range in which halos 
can be successfully detected. We treat this as an inde- 
pendent experiment from the previous section, beginning 
the comparison of ground- and space-based data afresh. 
In particular, we do not cut the Subaru data to the shal- 
lower depth discussed in §4.1[ to eliminate the last few 
systematic biases. The intent is not to align our two 
comparisons but rather to optimize each analysis as an 
independent experiment as would be the case if either 
were being undertaken as a self-contained survey. 

Unfortunately, even with the unprecedented invest- 
ment of HST time for the COSMOS survey, we can ex- 
pect the number of lensing-detected structures in this 
finite field to be modest. A t the Subaru depth, a surface 
density of ~5 halos deg -2 (|Mivazaki et al. 2007f) implies 
only around eight halos are likely to be found in the COS- 
MOS field. Thus we recognize in this comparison that 
the statistical significance of our results will be quite lim- 
ited. 

5.1. Residual Systematics 

First, we consider the B-mode signal. As discussed 
in the B-modes act as an independent realization of 
noise in the mass map, and locally highlight any prob- 
lems with the correction for PSF or other effects peculiar 
to the (two very different) instruments. Unsurprisingly, 
a visual inspection of Figure shows that the B-mode 
signal is significantly lower in our space-based data, with 
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fewer B-mode peaks. The overall noise level is reduced, 
and holes arising from masked foreground stars are also 
smaller and less frequent. In the ground-based maps, 
these create additional edges that lead to spurious effects 
during the Fourier transforms required by equation ([9]). 
The extended gaps are caused by difficulties modeling 
the PSF near the edge of the field of view, and could be 
eliminated in future surveys by more conservative tiling 
strategies. 

The southwest corner of the field has been troublesome 
throughout our analysis. This pointing was observed in 
slightly worse seeing, so the density of galaxies is reduced 
and the noise in the mass reconstruction is higher. 

5.2. Halo Detection 

The higher surface density of background galaxies from 
space also improves the reconstruction of the -E-mode 
"mass map" convergence field. The noise is lower and 
the angular resolution higher (although to aid compar- 
ison, both panels in figure [6] are smoothed to the same 
scale) . Several of the key features are qualitatively simi- 
lar but we are struck by the significant differences in the 
prominence of other mass peaks. To evaluate the robust- 
ness of detections, we shall now employ an automated 
peak-finding algorithm. 

Following (|Mivazaki et al. 20071 ). we smooth the con- 
vergence maps by a Gaussian kernel of rms width 1' and 
find local maxima with detection significance v > 4 (as- 
suming Gaussian errors on the shear measured within 
0.7' cells on the sky equal to the dispersion of those 
galaxy shears). Five peaks (marked A, B, C, D and 
K in figure [5]) are then identified in the ground-based 
data. However, two of these are near boundaries of the 
fiel d mask. Imposing the rigorous restrictions discussed 
by (|Mivazaki et al. 20071 ). we fin d that only peaks A , 
B and C remain (c.f. table 3 in iMivazaki et al. 20071) . 
All three are also detected in a space-based lens- 
ing analysis, the 3D distribution of galaxies and as 
extended sources in X -ray data (|Hasinger et al. 20071 
iFinoguenov et al. 2007|). Assuming the mass - luminosity 
relation adopted bv (jFinoguenov et al. 2007f l. the detec- 
tion threshold of this very deep X-ray data is well be- 
low that expected for lensing up to redshift ~ 1, so 
this acts as an ideal external arbiter (of course, X-ray 
mass-observable relations are somewhat uncertain) . The 
properties of the three clusters are summarized in ta- 
bleland demonstrate excellent agreement between the 
ground- and space-b ased data using the formalism of 
(TMivazaki et al. 20071 ). 

Cluster A (SJ J0959. 6+0231) is the most massive 
structure inside the COSMOS field, easily detected at 
many wavelengths. It appears to be in the process of 
a major merger, a nd has been studied individually by 
(jGuzzo et al. 20071) . who also obtained a spectroscopic 
redshift of z = 0.73. Cluster B (SL J1001.4+0159) is 
associated with an X-ray peak and overdensity of galax- 
ies at Zphot = 0.35. There is a second set of galaxies at 
Zphot = 0.85 within 2', which undoubtedly complicates 
the interpretation a little, but our results are consistent 
with this high redshift projection being a minor pertur- 
bation. Cluster C (NSC J100047+0 13912) is yet more 
local (z = 0.22. IMivazaki et al. 20071 ) and appears large 
on the sky. Only part of this cluster is inside the region 
of HST imaging, so the space-based signal is significantly 
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Fig. 6. — Convergence E-mode maps from the Subaru (left) and Hubble Space Telescope (right), af ter smoothing by a 1' Gaussian 
kernel. The data presented in the left panel are identical to that in figure 13 of (Miyazaki ct al. 2007), except that the field has been 
slightly truncated to match the right panel. Convergence is proportional to the total projected mass along a line of sight, modulated by 
the lensing sensitivity function (O plotted in figure[5] Contours are drawn at detection significances of 3<r, 4<r and 5<r, with dashed lines for 
underdensities. Clusters A, B, C and D are detected in both maps. O ther peaks E-L are on ly detected in one of the two. White enclosing 
circles denote clusters deemed "secure" by the rigorous standards of (Mivazaki ct al. 2007), and cyan circles denote "unsecure" clusters. 
The size of the circles shows the size of the smoothing kernel that maximises detection significance, enlarged by a factor of 2 for clarity. 




Fig. 7. — Convergence B-mode maps from Subaru (left) and Hubble Space Telescope (right) data. This is not produced by physical 
gravitational lensing, so deviations from zero include a combination of spurious effects from e.g. imperfect PSF correction, plus a realization 
of statistical noise. The smoothing scales, color ramps and contours are identical to those in figure [6] 
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weakened, and the mass is potentially underestimated by 
HST. 

To broaden our search, and test the limits of detectabil- 
ity, we additionally inve stigate the multi-scale procedure 
of (jHamana et al. 20 03). For this, we smooth the conver- 
gence maps with Gaussian filters of rms width between 
0.5' and 4', identifying local maxima inside the mask 
on each scale. For each peak with a detection signal to 
noise v > 4 on any scale, we record v and the smoothing 
scale that maximizes v. We also drop the restrictions on 
distance from the mask boundaries. This will increase 
the number of detected peaks, but at the expense of po- 
tentially introducing some spurious features. We then 
search for counterparts in the other data set, within 3' 
of detected peaks. 

With the above criteria, we identify four mass peaks 
common to both convergence maps (A, B, C and D). 
Cluster D is within 3' of an X-ray peak, and an over- 
density of galaxies at photometric redshifts z p hot = 0.93. 
This redshift is rather high for a lensing analysis, and 
it was flagged as "unsecure" by (|Mivazaki et al. 20071) 
because it is near a boundary in the image mask. How- 
ever, the tentative Subaru detection is strengthened by 
the confirmation from HST, and appears to be robust. 

Peaks E, F, G and H (also marked on Figure[6]) are seen 
only in the space-based map. The first three correspond 
to extended X-ray emission fr om clusters with masse s 
M 50 o between 2 - 4 x 10 13 M Q (|Finoguenov et al. 20071 ). 
Peak H is more massive (M500 = 1-8 x 10 14 M©), but is 
at very high redshift. All four of these peaks are real 
detections; however no counterparts within 3' are seen in 
the ground-based map, even down to v >3. Most likely, 
this is because of th eir lower mass and higher redshift 
(|Hamana et al. 20031 ). The detection of peak G was pre- 
vented in the ground-based data by a bright foreground 
star. 

Conversely, peaks I, J, K and L are detected only in the 
ground-based map. Peaks I and J are real, but lie just 
outside the HST imaging. There is an extended X-ray 
source at peak I, with unknown redshift, and a projection 
of two M 50Q « 2 x 10 13 M^ clusters at redshifts z = 0.40 
and z = 0.75 at peak J (|Finoguenov et al. 2007f ). In both 
cases, there is a weak, v < 3 signal in the HST data, 
fro m the wings of the clu ster. Peak K was detected in 
the (jMivazaki et al. 20071 ) analysis, but again flagged as 
"unsecure" because it is near the edge of a pointing. It 
does align with a slight, v < 3 detection in the space- 
based map, but there is no X-ray counterpart. This may 
be a spurious peak with chance coincidence, or perhaps 
a very distant object. Peak L appears to be spurious: 
such noise artifacts are more common near the edge of 
the field. 

In summary, to the extent that we can draw conclu- 
sions from such a small sample, there is very good agree- 
ment between the primary halo catalog drawn from the 
ground-based data and that independently found from 
the space-based data. Additional halos of lower mass 
and higher redshift are seen in the space-based catalog, 
and those located uniquely in the ground-based data can 
be understood in the context of either being outside the 
space-based region or close to its periphery. 

5.3. Halo Mass Estimation 



We now attempt to measure the total mass of each of 
the three halos (A-C) securely detected from both the 
ground and space. We assume that the clusters have an 
NFW density profile 

p(r) = 6 c p c /(r/r s )(l + r/r s ) 2 (10) 

(N avarro. Frenk fc White 19961 ). where 5 C is a function of 
the cluster's concentration c and scale size r s = 7"2oo/c, 
and r2oo is the radius within which the mean density 
is 200 times the critical densi ty. The shear profile o f 
an NFW cluster is derived by (jKing fc Schneider 2001h . 
We perform a maximum likelihood fit to the log(mass) 
and concentration parameters, using the shear measure- 
ments from all galaxies within 10' of the peak conver- 
gence signal, averaged in radial bins of 0.5'. It has 
been variously noted (J. Berge, S. Paulin-Hcnrikkson, 
private comm.) that fitting noisy data of individual clus- 
ters with an NFW profile does permit large (and there- 
fore massive) models with unnaturally low concentration 
values. To counter this effect, we impose a concentra- 
tion prior, using the lognormal distribution found for all 
hal oes in the Milleniu m Simulation as a function of mass 
by (jNeto et al. 20071 equations (5), (6) and figure 6). 
The resulting likelihood surfaces are shown in figure [5J 
with the effect of the prior being to close the bottom of 
the contours. Table [2] lists the best-fit masses and 68% 
confidence limits after marginalizing over concentration 
between 1 < c < 10. 

Although our common sample is small, there is an en- 
couraging agreement between the detailed properties of 
the clusters recovered from the ground and from space. 
For the higher redshift cluster A, our space-based data 
does put significantly tighter constraints on the mass 
and concentration than our ground-based data. How- 
ever, for the lower redshift clusters B and C, the results 
are satisfyingly similar. We note again that cluster C 
is partially outside the HST imaging. Since shears are 
only measured around one half of the cluster, the sta- 
tistical errors in the space-based analysis are larger and 
its mass could be underestimated. Certainly, for mas- 
sive clusters with redshifts 0.2 < z < 0.5, it appears that 
our ground-based depth and resolution is adequate. The 
main benefit of space-based imaging is in the measure- 
ment of lower mass halos and higher redshift clusters, 
plus the increased resolution to further investigate the 
distribution of their masses. 

6. DISCUSSION 

We have performed parallel weak lensing analyses of 
Subaru and Hubble Space Telescope imaging in the COS- 
MOS field. Our comparisons of the observed shear 
and convergence signals have revealed a number of is- 
sues, and suggest that such a study with real data use- 
fully complements the independe nt approach based on 
blind analyses of sim ulated data (jHevmans et al. 2 005a; 
iMassev et al. 2007cD . 

For statistical "cosmic shear" analyses, shear mea- 
surement with an existing ground based telescope, us- 
ing existing measurement techniques, can be achieved 
with less than 1% bias relative to higher resolution 
space-based data, for a galaxy surface density of 15 
axemin . One limitation of our approach is that we 
cannot check the performance of our space-based anal- 
ysis on the additional, small galaxies. At first sight 
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TABLE 2 

Cluster Masses 



Cluster RA Dec Redshift XMM mass HST mass Subaru mass 

10 14 M B 10 14 M o 10 14 M B 



A 149.917 2.515 0.73 1.90 ±0.05 23;" 13 
B 150.359 1.999 0.35 0.10 ±0.01 9±l 17 
C 150.184 1.657 0.22 1.01 ± 0.02 17±g 7 55. 




1 I ......I 

10 13 10 14 10 15 10 16 

Mass [M ] 

Fig. 8. — Best-fit mass and concentration index of three clusters 
in the COSMOS field, assuming NFW radial mass profiles. The 
contours show the 68%, 95% and 99% confidence regions obtained 
from Subaru data (dotted red) and Hubble Space Telescope data 
(solid black). 



the factor of ~3 shortfall in surface density seems in- 
consequential given the lower cost and improved areal 
mapping speed of existing ground-based cameras such 
as SuPrime-Cam. However, accompanying the brighter 
Subaru limit is a reduction in survey depth and hence 
the redshift distribution of background sources. More 
distant sources contain a larger signal, and a nar- 
rower range of redshifts also hinders tom ographical tests 
(jBacon et al. 20051 : iMassev et al. 2007bD . which tighten 
cosmological parameter constraints significantly. 

A key issue is whether this limiting depth is a 
fundamental one for all future ground-based cameras. 
PanSTARRS, VST, and even LSST each have signifi- 
cantly smaller primary mirrors than Subaru, so achieving 
even the S/N discussed here would require formidable ex- 
posure times. Most importantly, the deep infrared imag- 
ing that is required for photometric redshifts to enable 
tomographic analyses is likely to be difficult over large 
survey fields from the ground, because of increased sky 
background. Recent weak lensing analyses are limited at 
roughly the same level by uncertainty in galaxy shape 
measurement and photometric redshift estimation. 

Statistical measurements from the ground are also hin- 
dered by variable atmospheric seeing. Past experience 
has taught the authors that data collected in seeing worse 
than 0.8" is of little use for weak lensing analysis. The 
apparently rapid speed of data collection for our Subaru 
data belies the time spent waiting for better seeing, even 
with the excellent atmospheric conditions above Mauna 
Kea and the well-controlled dome seeing of Subaru. For 
this small-scale survey, we obtained exceptional quality 
imaging during a fortuitous observing window. The rele- 
vant quantity for larger-scale surveys in the future will be 
the time-averaged seeing quality, and the fraction of time 
spent with seeing better than 0.8". This is particularly 
true for surveys like Pan-STARRS and LSST, that plan 
to adopt a strategy of co-adding many shorter exposures. 
Their advantage is that the stacked images will achieve a 
near-uniform image quality, by virtue of the independent 
PSFs in each short exposure. This can then be tuned to 
the required image quality by rejecting a certain fraction 
of exposures. 

Variable seeing conditions is also of concern for the 
reconstruction of mass maps (c.f. Green et al. in prep.). 
Difficulties in the analysis of one pointing in the SW 
corner of the Subaru map result in a patchy recovery of 
large-scale structure; with more noise and a lower range 
of probed redshift in certain regions. 

Most importantly, however, the four most massive clus- 
ters out of eight detected from space are also detected 
from the ground - with one intriguing additional sig- 
nal and two more confirmed clusters just ouside the field 
of view observed from space. Reassuringly, the three 
clusters conservatively deemed "secure" by the indepen- 
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dent analysis of (|Mivazaki et al. 20071 ) have now been 
confirmed via space-based weak lensing and X-ray ob- 
servations. The physical properties of the four massive 
halos in common (A, B and C) are remarkably consis- 
tent whether derived from from ground- or space-based 
weak lensing. The measured masses and radial profiles 
of these clusters are consistent and, for the lower redshift 
clusters, the error bars are comparable. A Dark Energy 
Task Force "Stage 3" survey from the ground appears 
eminently feasible. 

A wide-field space-based platform would open up many 
new applications. Very important for statistical appli- 
cations is the increased redshift range of resolved back- 
ground galaxies. These not only contain a larger shear 
signal, but also more readily split into redshift bins for 
tomographic analysis. Three-dimensional analysis tech- 
niques will tighten constraints on cosmological parame- 
ters by factors of 3 — 5, and directly measure quantities 
that depend upon the properties of dark energy, like the 
growth of structure over cosmic time and the redshift- 
distance relation. They will also eliminate sources of er- 
ror due to the intrinsic correlations of galaxy shapes. Suf- 
ficiently good photometric redshifts require deep, wide- 
field near-IR imaging, and these are also realistically pos- 
sible over large surveys only from space. Recent weak 
lensing analyses with relativ ely shallow near-IR cover- 
age like ( jMassev et al. 200"7a| ) are limited to roughly the 
same degree by uncertainty in galaxy shape measurement 
and photometric redshift estimation. Full implementa- 
tions of cross-correlation cosmography will almost cer- 
tainly require deep near-IR imaging from space. Such 
advanced techniques will become particularly important 
as ground-based surveys expand to encompass the entire 
observable sky. 

The increased surface density of galaxies resolved from 
space also improves maps of the mass distribution. As 
shown in figures [6] and El the statistical noise and sys- 
tematic contamination in the -B-mode are significantly 
reduced. Eight clusters are detected in the COSMOS 
-E-mode signal without any contamination from spurious 
peaks. With the increased mass and spatial resolution of 
mass reconstructions from space, it becomes possible to 
detect halos the size of galaxy groups, as well as clusters 
over a wide range of redshifts - thus tracing their forma- 
tion, which is governed by the properties of dark matter 
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and the nature of gravity. Space-based data also crosses 
the threshold to mapping e ven filamentary large -scale 
structure in three dimensions (jMassev et al. 2007bl) . Ob- 
taining the detailed, 3D distribution of mass will be par- 
ticularly dmp_ortairtnearregions of interest like the Bullet 
cluster ([Clowe et al. 20061 ). where the small differences 
between the location of mass and baryons in a small 
patch of sky may yield the best possible information 
about the properties of dark matter. In this and other 
astrophysical phenomena, knowledge of the local mass 
environment and nearby large-scale structure is critical. 

Overall, we conclude that ground-based weak lensing 
surveys can perform several tasks remarkably well, with 
sufficiently small amount of systematic bias to easily jus- 
tify the next generation of dedicated ground-based sur- 
veys. Two dimensional statistical analyses will be able 
to produce order-of-magnitude improvements in weak 
lensing constraints, using proven hardware technology 
and software pipelines. On the other hand, a wide-field 
space-based imager would provide important control over 
some systematic effects, and open up many new applica- 
tions that are, at least currently, unachievable from the 
ground. For several of the most exciting techniques that 
will directly probe the nature of dark matter and dark 
energy, eventual space-based imaging is likely to be es- 
sential. 
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